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ABSTRACT 

We present a Monte Carlo calculation of the astrophysical rate of the 
15 0(a, 7) 19 Ne reaction based on an evaluation of published experimental data. 
By considering the likelihood distributions of individual resonance parameters 
derived from measurements, estimates of upper and lower limits on the reac- 
tion rate at the 99.73% confidence level are derived in addition to the recom- 
mended, median value. These three reaction rates are used as input for three 
separate calculations of Type I x-ray bursts using spherically symmetric, hydro- 
dynamic simulations of an accreting neutron star. In this way the influence of 
the 15 0(a, 7) 19 Ne reaction rate on the peak luminosity, recurrence time, and as- 
sociated nucleosynthesis in models of Type I x-ray bursts is studied. Contrary to 
previous findings, no substantial effect on any of these quantities is observed in 
a sequence of four bursts when varying the reaction rate between its lower and 
upper limits. Rather, the differences in these quantities are comparable to the 
burst-to-burst variations with a fixed reaction rate, indicating that uncertainties 
in the 15 0(o;, 7) 19 Ne reaction rate do not strongly affect the predictions of this 
Type I x-ray burst model. 

Subject headings: accreting neutron stars, Type I x-ray bursts, nuclear reaction rate 
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Introduction 



Type I x-ray bursts are thought to arise from thermonuclear runaways occurring in a 



thin layer on the surfaces of neutron stars 



companion stars in close binary systems (jWoosley fc Taam 



accreting hydrogen- an c 



1976 



helium-rich matter from 



Joss 



19771 ). Accreting 



neutron stars in binary systems stably fuse hydrogen into helium via the temperature- 
independent hot CNO cycles at high accretion rates. At low accretion rates they exhibit 



x-ray bursts due to the explosive burning of H or He o 
which approximately 10 39-40 erg is liberated (jCumming 



typic al duration 10-100 s during 



20041 ). If the accreted material is 
rich in H, the rp proce ss occurs after the burst ignite s, potentially leading to the synthesis o f 



nuclei up to A « 100 ( ISchatz et al. 



20011 ) or beyond (IKoike et al 



2004 : 



Elomaa et al 



2009|). 



The light nuclei that participate in the CNO cycles and the heavier nuc 



are linked by the breakout reactions 15 0(a,7) 19 Ne and 18 Ne(a,p) 21 Na ([Wallace fc Woosley 



ei of the rp process 



1981 



Wiescher et al. 



19991 ) . The latter reaction becomes important at higher temperatures 
than does the former; hence 15 0(a,7) 19 Ne has a larger effect on the hydrogen mass fraction 
at the time of burst ignition than does 18 Ne(a,p) 21 Na, thereby more strongly influencing 
the dynamics of the burst. For this reason, t he effect of this rate on x-ray burst models 



has been the subj ect of several investigations ( IFisker et al. 



Fisker et al. 



20071 ). 



20061 : 



Cooper fc Narayan 



20061 : 



While the 15 0(a,7) 19 Ne reaction rate has not been measured directly, intense and 
steady experimental efforts have enabled measurements of the decay properties of the states 
in 19 Ne relevant to the reaction rate at x-ray burst temperatures. These states are just 
above the a threshold which lies at an excitation energy of 3.53 MeV. As they lie below the 
proton and neutron thresholds, the states of interest decay only by a or 7 ray emission. 
Hence the resonant rate of 15 0(a, 7) 19 Ne can be calculated from the radiative and a widths, 
or equivalently the lifetimes and a-decay branching ratios of the relevant states. 
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In this paper we assess the experimental data on the decay properties of states just 
above the a emission threshold in 19 Ne and use a Monte Carlo technique to evaluate the 
reaction rate. With the results of this analysis upper and lower limits on the reaction rate 
at the 99.73% confidence level are derived in addition to the recommended, median value. 
These three reaction rates are used as input for three separate calculations of Type I x-ray 
bursts using spherically symmetric, hydrodynamic simulations of an accreting neutron 
star. In this way the influence of the 15 0(o;, 7) 19 Ne reaction rate on the peak luminosity, 
recurrence time, and associated nucleosynthesis of Type I x-ray bursts is studied and 
compared with previous findings. 

2. Decay Properties of 19 Ne States and the Resonant Rate of ls O(a, 7) 19 Ne 

The thermonuclear rate of the 15 0(o;, 7) 19 Ne reaction at temperatures below 2 GK, 
the relevant range for Type I x-ray bursts, is dominated by the contributions of resonances 
corresponding to states in 19 Ne with excitation energies between 4 and 5 MeV. These states 
lie above the a decay threshold but below the thresholds for neutron or proton decay. 
Hence they decay only by a and 7 emission and the total decay width T is the sum of 
these partial decay widths, i.e., r = T a + T 7 . As these resonances are narrow compared to 
the spacing between adjacent levels with the same spin and parity, the contributions of the 
different resonances add incoherently and can be measured separately. Even in the absence 
of a direct measurement, the strength of a resonance can be inferred from measurements of 
the radiative and a widths of a state. Equivalently, measurements of the mean lifetime of a 
state r and its a-decay branching ratio B a = r Q /T suffice to determine the strength of a 
resonance. The resonance strength W7 of a 19 Ne state with angular momentum J formed 
by the capture of an a particle by 15 O is given by 



Co>7 = 



2j+i r Q r. 



2J + 1 (B a F)(l - B a )T _ 2J + 1 



B a (i-B a )r. (1) 



(2j a + i)(2ji 50 + i) r 



2 r 2 
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Taking account of the fact that the mean lifetime and total decay width are related by 
T = h/r, this implies that 

(2J+l)B a (l-B a )H 
u l = — • ( 2 ) 



2r 



over the last twenty years ( 


Maenus e 


t al. 


1990 




Laird et al. 


2002; 


Davids et al. 


2003a 




Davids et al. 


2003; 


Rehm et al. 


2003 




Visser et al. 


2004; 


Tan et al. 


2007; 


Tan et al. 


2009) 



The results of a numb er of B n measurements of states in 19 Ne have been publish ed 



In addition, the lifetimes of these states have been measured 



over the past five years (ITan et al. 



2005 



Kanungo et al. 



2006; 



jy two independen t groups 



Mvthili et al. 



20081 ). Here we 



critically evaluate the published measurements. 



2.1. Mean Lifetime Data 

All published lifetime measurements for states with excitation energies up to 4.602 
MeV in 19 Ne are shown in Table [D At least two measurements of the lifetimes of all the 
known states with excitation energies between 1536 keV and 4602 keV have been published 
and agree. This agreement gives us a measure of assurance that the values are reliable, 
though apart from the state at 4.03 MeV, the only Doppler shifted 7-ray sp ectra for state s 



above the a threshold that have been published in the literature appear in 



Mvthili et al 



(120081 ). For the states lying above 4602 keV, no published lifetime data exist, likely because 
the states decay predominantly by a decay and therefore were not observed in Doppler 
shift attenuation measurements. We use analog state information to deduce information 
on the radiative widths of the 4.71 and 5.09 MeV states. Although there are known 



states above 5.092 MeV and below t 



ratios consistent with 1 (jRehm et al 



re pro ton threshold, they all have a decay branching 



20031 ) and their high excitation energies render them 
unimportant contributors to the thermally averaged reaction rate at Type I X-ray burst 
temperatures. Hence we make no attempt to estimate the lifetimes of these states. 
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2.2. Spin and Parity Assignments 

In the relevant excitation energy range, the spins and parities of most of the states are 
well known. The exceptions are two of the negative parity states, the 4144 keV state with a 
tentatively as s igned spin of 9/2 and the 4200 keV state tentatively assigned to be spin 7/2. 



Garrett et al. 



( 197: 



2) initi 



state 9/2 state in t he context of 
on their DWBA analyses iGarrett et al. 



tially discussed the 4144 keV state as a 7/2 state and the 4200 keV 
;he known rotational bands in 19 F and 19 Ne. Based 



( 119721 ) suggested swappi ng the spin assignm ents of 



these two states, a suggestion tentatively adopted by evaluators ( ITilley et al 



1995|). 



We compare the spins of the 19 Ne states with the spins of the presumed isobaric analog 
states in 19 F by examining the reduced transition probabilities and the branching ratios 
of putative analog states. The isoscalar component dominates strong E2 transitions. The 
4033 keV state in 19 F undergoes a strong E2 transition. Depending on the spin, either the 
4144 keV state in 19 Ne or the 4200 keV state under goes an E2 transition. Both the 3999 
and 4033 keV states in 19 F decay to the 5/2~ state at 1346 ke V. We calculate the r educed 



transition probabiliti es of the 19 Ne state s usi ng the lifet i mes o f 



transition energies of 



Tilley et al. 



fll995f ) and 



Tan et al. 



Mvthili et al. 



()2008|) and 



((20051) . In 19 F, the 3999 keV state 



decays to the 1346 keV state via an Ml transition while the 4033 keV state undergoes an 
E2 transition. These transition strengths are tabulated in the first row of Table El The 
4144 keV and 4200 keV states in 19 Ne decay to the 5/2" state at 1508 keV. If we accept the 
tentative spin assignments of 9/2~ and 7/2~ respectively, the 4200 keV state undergoes an 
Ml transition and the 4144 keV state an E2 transition. The resulting reduced transition 
probabilities appear in the second row of Table El The reduced transition probabilities 
calculated with the alternate spin assignments are shown in the third row of Table El where 
we see improved agreement with the 19 F B(E2) value. 



Although the reduced transition probabilities suggest that perhaps the spins should be 
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swapped, th e mea sured 7 decay branching ratios support the tentative spin assignments of 



Tillev et al. 



(119951 ) for the 4144 keV and 4200 keV states in 19 Ne. The 9/2" state at 4033 
keV in 19 F decays exclusively to the 1346 keV state; the 4144 keV state in 19 Ne decays 
exclusively to the 1508 keV, suggesting that the tentative spin assignment is justified. The 
4200 keV state in 19 Ne and the 3999 keV state in 19 F have comparable branching ratios of 
80% and 70% to the 1508 keV state in 19 Ne and the 1346 keV state in 19 F respectively. It 
is necessary that the spins of the 4144 and 4200 keV states be measured experimentally in 
order to precisely constrain their contributions to the reaction rate, th ough this is not the 



largest uncertainty; for now we adopt the tentative spin assignments of 



Tillev et al 



(1995) 



2.3. a Decay Branching Ratio Data 



Independent measurements of the a decay branching ratio have yielded good agreement 
for the high lying states, but measurements of the low lying states have been controversial. 
All published measurements of whic h we are aware are sh own in Table [3] and Table HI The 



first measurements were reported by 



Magnus et al. 



(119901 ) and are shown in Table [3j The 



lowest energy state whose a decay branching ratio was reported in this measurement was 
the 4.38 MeV state, for which B n = 0.044 ± 0.032 wa s given. Examination of the timing 



spectrum for this state shown in 



Magnus et al. 



(Il990l ) reveals no peak and leads us to 



conclude that, as one would expect a priori, a 1.4a excess over background such as this does 
not warrant a detection claim; these data can only be used to constr a in _E? 7 from above. 



In contrast, all the other B a measurements reported in 



Magnus et al. 



(119901 ) appear to be 



based on statistically significant excesses observed over well-determined backgrounds. 



Laird et al. (2002 



that of 



Magnus et al 



mea sured B a for a single state and found a value consistent with 



(|l990h 
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The measurement reported in iDavids et al.l (120031 ) was the most sensitive made to date 
of the a decay of the 4.03 MeV state owing to the small, well-measured background and the 
selectivity of the p( 21 Ne,£) 19 Ne reaction used to populate it. Contrary to the claim made in 



Tan et al. 



(120091 ). net a decay events above the estima ted background wer e observed at the 



Davids et al 



(120031 ). However this 



location of the 4.03 MeV level in the measurement of 
small excess was not deemed statistically significant and the measurement was reported as 
an upper limit. No excess above the background was observed at the location of the 4.38 
MeV state and therefore an upper limit on B a was also reported for this state. We note that 
this 90% confidence level upper limit on B a for the 4.38 MeV state is more than a factor 
of ten smaller than the ce ntral value of the sta tistically insignificant detection of a decay 



reported from this state in 



Magnus et al. 



(119901 ) . illustrating the peril of relying on gaussian 



representations of the likeli hood distribut i ons o f such measurements. The B n values for 



the ot her states reported in IDavids et al. 



(1990) 



( 120031 ) are consistent with those of iMagnus et al 



Rehm et al. 



(120031 ) report B a measurements of three states in 19 Ne but the excitation 
energy resolution was around 220 keV and was insufficient to resolve adjacent levels, 
rendering the fitting used to deduce the two B a central values and associated uncertainties 
less than wholly satisfying. Nevertheless there was no background around the 4.03 MeV 
level, permitting a firm upper limit on B a for this state to be established. 



Visser et al. 



( 120041 ) 



Problems with the electronics used in the measurement of 
rendered the normalization of the B a values reported there questionable. In particular, 
the measurements were normalized to the coincidence efficiency determined based on B a 
measurements of the 5351 keV level. The measurement yielded different values of B a for 
this state, all of which differed substantially from the accepted value of 1, depending on the 
different electronics gating configurations used in the measurement. Out of concern that 
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the normalization is not completely understood, we have avoided using these data. 



Although the relev ant states were po pulated and separated with good energy resolution 



in the measurement of 



TanetaL 



(120071 ) . the large background was not persuasively 



modelled for the states below 4.55 MeV from which a decay was putatively detected. Since 



the backgrounds were not well determined in the 4.03, 4.1 



decay energy spectra shown in Figure 3 of 



Tan et al 



/4.20, and 4.38 MeV state a 



(120071 ). the excesses claimed to be 



observed above them are not quantitatively reliable. Inspection of the timing spectra for 
thes e putative a deca y detections from states below 4.55 MeV in 19 Ne, shown in Figure 



6 of 



Tan et al. 



(120091 ). reveals no convincing evidence of an a decay peak from the 4.03, 



4.14, 4.20, or 4.38 MeV states. In contrast the B a values for states lying at and above 4.55 
MeV appear reliable as far as the background detern 
concerned, and they agree with the measurements of 



MeV appear reliable as far as the background deter mination and s i gnal t o noise ratio are 

~x3). 



Magnus et al. 



(I199C 



3. Monte Carlo Reaction Rate Calculation 

We report here the results of a million event Monte Carlo simulation of the thermally 
averaged rate of the 15 0(a,7) 19 Ne per particle pair as a function of temperature. 
This reaction rate is the incoh erent sum of the nonresonant rate, as calculated by 



Dufour fc Descouvemont 



(|2000| ). and the contributions due to the 19 Ne states at 4.03, 4.14, 
4.20, 4.38, 4.55, 4.60, 4.71, and 5.09 MeV. The nonresonant contribution is negligible above 
0.2 GK. While the nonresonant rate is fixed at a central value due to its totally insignificant 
impact on the error budget, the resonant contributions are calculated by drawing randomly 
from the likelihood distributions of the resonance parameters. The energies and spins of 
the levels adopted here are shown in Table [5J The contribution due to each resonance was 
calculated as f ollows . The e nergy and la uncertainty of each state was taken from the 



compilation of 



Tilley et al. 



(119951 ) or calculated by using a weighted average of this and 
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the value from iTan et al.l ( 12005 



discrepant data as described in 



, including a 



Cyburt et al. 



discrep ancy error to account for potentially 



(1200 ll ). The likelihood distributions for the 



level energies are taken to be gaussians characterized by the means and standard deviations 
given in Table [5j 



The lifetime of the 4.03 MeV state is drawn from the likelihood distributions 
quantity reported by 



Tan et al. 



(120051 ) 



Kanungo et al. 



O2006|), and 



Mvthili et al. 



or this 



( I2008h 



with the weight of each measurement assigned according to the inverse of its variance. The 
resulting likelihood distribution realized in the one million event Monte Carlo simulation is 
shown in Figure [H For B n of this st ate we use the likelihood distribution derived from the 



measurement of 



Davids et al. 



(120031 ) , shown in Figure [2j 



Prior to the measurement of 



Davids et al. 



for the 



(120031 ). the best information on 
4.03 MeV state came from measurements of a transfer to the 19 F analog state (iMao et al. 



19951 ). The uncertainties in using isospin symmetry to deduce reduced alpha widths o f 



state s near threshold from their analog states are fairly well quantified (IFortune et al. 
20101 ). but systematic uncertainties due to the reaction mechanism remain. Fortunately 
there are reliable experimental data on B a for the 4.03 MeV state and we need not rely 
on the analog state. Moreover the most likely value of r a for the 4.03 MeV state derived 
from the Monte Carlo si mulation is 8 MeV , perfectly consistent with the value of 8.8(14) 



/xeV given (for q = 9) in 



Mao et al. 



(1995) 



The mean lifetimes of the 4.14 and 4.20 



of 



Tan et al. 



(120051 ) and 



Mvthili et al 



MeV states are drawn from the measurements 
2008f ) as described above. As we are not convinced 
of the reliability of the single B a measurement for these two states reported in the literature, 
we calculate the a widths of these states according to 

3h 2 



(3) 



where \x is the reduced mass of the 15 0+a system, a the channel radius, P the Coulomb 





5 10 15 20 25 30 

Mean Lifetime of 4.03 MeV State (fs) 

Fig. 1. — Likelihood distribution of the mean lifetime of the 4.03 MeV state in 19 Ne based 
on the measurements of Tan et al. 2005, Kanungo et al. 2006, and Mythili et al. 2008. The 
weight assigned to each measurement is proportional to the inverse of its variance. 
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Table 1: 19 Ne level mean lifetime measurements and la uncertainties in fs. Statistical and 

systematic errors are listed separately in the last column. 

Level Energy (MeV) Tan et al. (2005) Kanungo et al. ( 20061 Mvthili et al. (2008) 



1.54 16 ±4 19.1^ ±1.1 
4.03 13tl n -3 6.9 ±1.5 ±0.7 
4.14 18±§ 14.01^0 ±1-2 
4.20 43^ 2 381^ ± 2 
4.38 5±| 2.9 ±1.4 ±0.6 

4.55 1511 1 18.7±i;g ±2.2 
4.60 7±\ 7.6±|;5 ± 0.9 




0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 
4.03 MeV State a Decay Branching Ratio 



Fig. 2. — Likelihood distribution of the a decay branching ratio of the 4.03 MeV state in 
19 Ne based on the measurement of Davids et al. 2003. 
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Nucleus J 71 " Energy Level 5 (Ml) Nucleus J n Energy Level B(E2) 
(keV) (MeVfm 3 ) (keV) (MeV fm 5 ) 

19 F |- 3999 0.00171^005 19p |~ 4033 90±2 ° 



19 Ne (|)- 4200 0.0008±g;ggg| 19 Ne (§)" 4144 460^° 



19 Ne (I)- 4144 0.0024t°;°oo9 19 Ne (§)" 4200 150±gg 



Table 2: Reduced Ml and E2 transition probabilities for two isobaric analog states in 19 F 
and 19 Ne. The first row shows the transition strengths for two states in 19 F. The second 
row contains the 1 9 Ne transition strengths with the tentatively adopted spin assignments of 



Tilley et al. 



( 119951 ) for the 19 Ne analog states and the third row gives the transition strengths 
with the alternate assignment of the 19 Ne spins. 



penetrability, and Q\ the reduced a width (ITeichmann fc Wiener 



efforts to determine Q 2 a for the 19 F analog states (jde Oliveira et al. 



1952). Despite some 



1996 



Mao et al. 



1996|), 



in our judgement no reliable experimental determinations of these reduced widths are 
available in the literature, so we draw them from distributions of our own construction. 
We conjecture that these reduced widths can be adequately represented by a log normal 
distribution having a mode of Q\ = 0.005 and standard deviation of its logarithm of 1.5. 
In a survey of ff 1 valu es of unbound states in the a nuclei 24 Mg, 28 Si, 32 S, 36 Ar, and 40 Ca, 



Longland et ai- 



de Oliveira et al. 



( 2010h 



find a mean value of 0.01. For comparison, the Q 2 a value given in 



(1996) (without any uncertainty) is 0.135 for both states. We chose a 
channel radius of a = 5.8 fm, corresponding to r = 1.45 fm; we investigated different 
channel radii of 5.6 and 6.0 fm and found that the resulting penetrabilities changed by 
~ 50%, which is totally insignificant compared with the uncertainty due to the reduced 
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Table 3: 19 Ne level a decay branching ratio B a measurements and la uncertainties. Upper 



limits are cited at the 90% confidence level. 


Level Energy (MeV) 


Maeniis et al. (1990) 


Laird et al. (2002) 


Davids et al. (2003) 


4.03 






< 4.3 x 10~ 4 


4.38 


0.044 ±0.032 




< 3.9 x 10~ 3 


4.55 


0.07 ±0.03 




0.16 ±0.04 


4.60 


0.25 ±0.04 


0.32 ±0.03 


0.32 ±0.04 


4.71 


0.82 ±0.15 




0.85 ±0.04 


5.09 


0.90 ±0.09 




0.90 ±0.06 



Table 4: Ne level a decay branching ratio B a measurements and la uncertainties. Upper 
limits are cited at the 90% confidence level. 



Level Energy (MeV) Rehm et al. (2003) Visser et al. (2004) Tan et al. (2009) 



4.03 


< 6 x 10~ 4 




2.9 ±2.1 x 10~ 4 


4.14 and 4.20 






1.2 ±0.5 x 10~ 3 


4.38 


0.016 ±0.005 




1.2 ±0.3 x lO" 3 


4.55 




0.06 ±0.04 


0.07 ±0.02 


4.60 




0.208 ±0.026 


0.26 ±0.03 


4.71 






0.80 ±0.15 


5.09 


0.8 ±0.1 


0-75±g5? 


0.87 ±0.03 
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width. The likelihood distribution for the B a value of the 4.20 MeV state resulting from 
this calculation is shown in Figure [3J 

We treat the 4.38 MeV state in the same way as we do the 4.03 MeV state, with 



the exception being that only two independent measurements of r exist. The 



distribution of B Q for this state is taken from the measurement of 
is shown in Figure HI 



Davids et al. 



telihood 



(|2003j) and 



of 



The mean lifetimes of the 4.55 and 4.60 MeV states are drawn from the measurements 

values of the 4.55 and 4.60 MeV states 



Tan et al. 



(120051 ) and 



Mvthili et al 



are drawn from the measurements of 



(2008 ). TheB 



fl2003h . and 



Tan et al. 



Magnus et al. 



(1990) 



Laird et al. 



( 20021 



Davids et al. 



( 120091 ). with the weight of each measurement proportional to the 



inverse of its variance as usual. 

In the cases of the 4.71 and 5.09 MeV states for which lifetime data are lacking we 
assume that isospin symmetry is valid and equate the reduced transition probabilities of 
these states with t heir analog states in 19 F. Utilizing the lifetime measurements compiled i n 



Tilley et al. 



(119951 ). the 7 decay branching r atio measurement o 



and the resonance strength measurement of 



Wilmes et al. 



Pringle fc Vermeerl (119891 ) 



widths of the 4.71 MeV and 5.09 MeV states are 46± " meV and l in+] 10 



'he B„ values are drawn from the measurements of 



(I2003h . and 



Tan et al 



(I2009[ ). 



(120021 ). we find that the radiative 
0l«n° m eV respective ly. 



Magnus et al. 



(1990) 



Davids et al 



Using the resonance parameter likelihood distributions described above, we calculated 
the thermally averag ed rate of the 15 0(a,7) 19 Ne per particle pair as a function of 
temperature one million times at a range of temperatures between 0.1 and 2 GK. Analysing 
the results, we determine upper and lower limits on the reaction rate at the 99.73% 
confidence level and adopt the median of the calculated reaction rate distribution as our 
recommended value. The calculated reaction rates are shown in Figure [5] and numerical 
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C/3 

a 

o 

1 

•T-H 

00 



120000 



100000 




0.000 



0.002 0.004 0.006 0.008 

4.20 MeV State a Decay Branching Ratio 



0.010 



Fig. 3. — Likelihood distribution of the a decay branching ratio of the 4.20 MeV state in 
19 Ne. We compute this using R-matrix formulae assuming the reduced a width is given by 
a log normal distribution peaked at 0.005 and having a standard deviation of its logarithm 
of 1.5. 
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Table 5: 19 Ne level angular momenta and energies with la uncertainties adopted in this work. 



Level Energy (keV) 


j w 


Reference (s) 




4034.3(9) 


3/2+ 


Tillev et al. 1995): Tan et al. 


(2005) 


4143.4(8) 


(9/2)- 


Tillev et al. (1995): Tan et al. 


(2005) 


4199.7(16) 


(7/2)" 


Tillev et al. (1995): Tan et al. 


(2005) 


4377.9(7) 


7/2+ 


Tillev et al. 1995): Tan et al. 


C2005) 


4547.8(10) 


3/2- 


Tillev et al. 1995); Tan et al. 


f2005) 


4601.7(9) 


5/2+ 


Tillev et al. 1995): Tan et al. 


f2005) 


4712(10) 


5/2- 


Tillev et al. 1995) 




5092(6) 


5/2+ 


Tillev et al. 1995) 






Fig. 4. — Likelihood distribution of the a decay branching ratio of the 4.38 MeV state in 
19 Ne based on the measurement of Davids et al. 2003. 
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values are given in Table |6j 



m 



Graphical compar i sons o f the Monte Carlo r eactio n rates to the reaction rates estimated 



Caughlan k Fowleri (119881 ) and 



Fisker et al. 



(120061 ) are shown in Figure [6] and Figure 



[7] respectively. In the temperature range from 0.1-2 GK the Monte Carlo low e r lim it on 



the reaction rate ranges between 2% and 91% of the lower limit of 
Monte Carlo central value lies between 0.71 and 1.04 times that of 
the Monte Carlo upper limit is between 3.2 and 9.4 times that of 



Tan et al. 



Tan et al. 



Tan et a 



Tan et al. 



(2009), the 
(12009k and 
(12009L Hence 



(I2OO9J) fall well 



the central value and limits on the reaction rate estimated by 
within the 99.73% CL limits on the reaction rate obtained from the Monte Carlo simulation. 
Comparing to the recent reaction rate estimates of 



Iliadis et al. 



( 120 lOl ). we find that over 



the same temperature range of 0.1- 2 GK our Monte C arlo lower limit on the reaction rate is 



between 0.03 and 3.6 times that of 
3.3 times that of 



Iliadis et al 



iadis et al. 



liadis et al 



( 120101 ). our central value is between 1.7 and 



( 120101 ). and our upper limit is between 7 and 14 times that of 



( 120101 ). However, throughout the entire temperature range characterizing the 



initial phase of the thermonuclear runa way (below 0. 7 GK 



the reaction rate is smaller than that of 



Iliadis et al. 



(I2010f ). 



our Monte Carlo lower limit on 



Given the large uncertainties in the B a values of the low lying states, it is impossible 
to precisely specify the fraction of the total reaction rate accounted for by the resonant 
contribution of any particular state. But in the Monte Carlo approach taken here it is 
possible to find the likelihood of a given state accounting for any definite fraction of the 
total reaction rate. Figure E] shows the likelihood distributions of the fraction of the total 
reaction rate at 0.5 GK accounted for by the resonant contributions of the 4.03 and 4.14 
MeV states. It is possible to say that at this temperature and indeed below 0.6 GK, the 
4.03 MeV state likely dominates the reaction rate. 
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Temperature (GK) 



Fig. 5. — Product of the Avogadro constant and the thermally averaged rate of the 
15 0(a, 7) 19 Ne per particle pair as a function of temperature. Shown are the median and 
99.73% confidence level upper and lower limits of the probability distribution from the Monte 



Carlo simulation. 
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Table 6: Product of 


the Avogadro constant 


and the thermally averaged rate of 


15 0(a, 7) 19 Ne reaction per particle pair in cm 3 


mol 1 s 1 




Temperature (GK) 


99.73% CL Lower Limit 


Median 


99.73% CL Upper Limit 


0.1 


1.1 x 10~ 25 


1.0 x 10~ 24 


7.7 x 10" 24 


0.2 


6.1 x 10~ 15 


1.8 x 10~ 12 


1.4 x lO" 11 


0.25 


2.7 x 10~ 12 


4.6 x 10~ 10 


3.6 x 10~ 9 


0.3 


1.6 x 10- 10 


1.8 x 10~ 8 


1.4 x 10~ 7 


0.35 


3.3 x 10~ 9 


2.3 x 10~ 7 


1.8 x 10~ 6 


0.4 


3.2 x 10~ 8 


1.6 x 10~ 6 


1.2 x 10~ 5 


0.45 


2.1 x 10~ 7 


7.0 x 10~ 6 


5.1 x 10^ 5 


0.5 


1.0 x 10~ 6 


2.3 x 10~ 5 


0.00017 


0.6 


1.3 x 10~ 5 


0.00015 


0.0012 


0.7 


0.00011 


0.00061 


0.0069 


0.8 


0.00066 


0.0021 


0.033 


0.9 


0.0028 


0.0066 


0.11 


1.0 


0.0088 


0.018 


0.31 


1.25 


0.074 


0.13 


2.0 


1.5 


0.31 


0.53 


7.4 


1.75 


0.85 


1.5 


18 


2.0 


1.78 


3.1 


35 
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10 =- 



1 r- 



0.1 - 



0.01 



0.001 




0.2 0.4 0.6 



0.8 1 1.2 1.4 
Temperature (GK) 



1.6 1.8 



Fig. 6. — Ratio o f the Monte Carlo reaction rate to the reaction rate estimated in 

re solid line is the ratio of the recommended reaction rate 



Cauehlan & Fowlei 


f 


1988 


). 1 


to that of 


Cauehlan & Fowler 



(119881 ) and the shaded region is allowed at the 99.73% confi- 



dence level. 
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10 =- 



1 r- 



0.1 - 



0.01 £ 



0.001 




0.2 0.4 0.6 



0.8 1 1.2 1.4 
Temperature (GK) 



1.6 1.8 



Fig. 7 . — Ratio of the Monte Carlo reaction rate to the reaction rate estimated in 



(2006). The solid line is the ratio of the recommended reaction rate to that of 



Fisker et al. 



FiskeretaL 



(120061 ) and the shaded region is allowed at the 99.73% confidence level. 
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4.03 MeV State 
4.14 MeV State 
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0.2 0.4 0.6 0.8 

Fraction of Total Reaction Rate at 0.5 GK 



Fig. 8. — Likelihood distributions of the fraction of the total reaction rate at 0.5 GK ac- 
counted for by the 4.03 and 4.14 MeV states. Below 0.6 GK, the reaction rate is likely 
dominated by the contribution of the 4.03 MeV state with small contributions from the 4.14 
and 4.20 MeV states. 
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4. Simulations of Type I X-Ray Bursts 



In order to assess the impact of the current uncertainties in the 15 0(a, 7) rate on Type 
I x-ray burst (XRB) models, we have performed a new set of simulations with SHIVA, a 
one dimensional (spherically symmetric), hydrodynamic, implicit, Lagra ngian code used 



modeling 


: stell 


(Jose et al. 


2010) 



bursts, with identical input physics except for the specific choice of the 15 0(a, 7) 19 Ne rate 
have been computed: the lower limit reaction rate dubbed Model A, the recommended rate 
Model B, and the upper limit Model C. In all cases the accretion of solar-like matter with 
hydrogen mass fraction X = 0.7048, helium mass fraction Y = 0.2752, and metallicity Z = 
0.02 onto a 1.4 M & neutron star (with initial luminosity Li = 1.6 x 10 34 erg s -1 = 4.14 Lq), 
at a mass accretion rate M acc = 1.75 x 10~ 9 M & yr _1 (corresponding to 0.08 MEdd) has 
been adopted. The Eddington accretion rate MEdd is the accretion rate that produces a 
luminosity at which the force due to radiation pressure equals the local gravitational force. 
All metals are initially assumed to be in the form of 14 N, which is justified considering 
the rapidity with which the CNO nuclei are converted to 14 N at the high temperatures 
reached early in the burst and the negligible energy this releases compared with accretion 



( jWoosley et a 



computed by 



2004) . The m odel 



Woosley et al. 



adopted in this work is qualitatively similar to model zM 



( 120041 ) in the framework of the ID, 



code KEPLE R, as well as to 



that whereas 



Woosley et al 



;he m odel computed by 



i vdro d ynamic, implicit 



Fisker et al. 



(120061 ) . We note however 



( 120041 ) assume a value of 10 km for the (Newtonian) neutron 
bllowing the integr ation of the neutron star 



Fisker et al. 



(120061 ). in a general relativistic 



star radius, our model yields a value of 13.1 km 
structure. In turn, the calculations reported by 
framework, relied on a radius of 11 km for the same neutron star mass. Differences in 
the neutron star size and thereby surface gravity may affect the strength of the explosion 



through the accreted mass, peak temperature, nucleosynthesis, etc. to some extent. 
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The code has been linked to a fully updated nuclear reaction network containing 324 
nuclides from 1 H to 107 Te and 1392 nuclear processes, and includes the most important 



charged particle induced reactions occurring am ong the nuclides bet ween and 



107' 



Parikh et al. 



fl2008|) and 



Jose et al. 



e as 



(boioh 



well as their corresponding inverse processes (see 
for details). Neutron captures are not taken into account since they play a very minor role 
in XRB nucleosynthesis. SHIVA uses a time-dependent formalism for convective transport 
whenever the characteristic convective timescale becomes larger than the integration time 
step. Par tial mixing between adjacent convective shells is treated by means of a diffusion 



equation (IPrialnik et al. 



19791 ) . Accre tion is computed 



a constant number of envelope shells (IKutter fc Sparks 



jy red istributing material through 



19801 ): a thin envelope, containing 



1.1 x 10 18 g of material (less than 1/1000 of the total envelope mass accreted during 
the first bursting episode), equally distributed through all the envelope shells, is chosen 
as an initial condition. The model is then relaxed using a few very large time steps to 
gu arantee hy d rosta tic equilibrium. More details on the adopted input physics can be found 



in 



Jose et al. 



feoioh . 



5. Results and Discussion 

A sequence of four consecutive bursts has been computed for all three models reported 
in this work. Summaries of the gross properties of the different sequences are given in 
Figure [9j Table [7J Table [8] and Table [9j For each model, the mean mass fractions of all 
nuclides i that are stable or have a half-life >1 h and are present with a mass fraction 
Xi > 10~ 9 in the envelope at the end of the fourth burst are given in Table [101 

As a framework for our discussion, we present first a summary of the results for Model 
B, in which the recommended 15 0(a,7) 19 Ne rate has been adopted. As shown in Tabled 
peak temperatures T pea }. ~ 1.1 — 1.3 GK and peak luminosities L pea k ~ (1 — 2) x 10 5 L & 
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are reached in this model. Recurrence times between bursts around r rec ~ 4.5-6 hr and 
a ~ 28 — 50, have also been obtained. The ratio between persistent and burst luminosities, 



a, is defined as 



jt+T rec ^ df 



a 



/;™L(t) dt 



(4) 



"t'+TQ.Ql 
t> 

with the latter integrated over the time during which the burst exceeds 1% of 

its peak luminosity, To.oi- During the interburst period, the accretion luminosity, 

L aC c = GMM/R ~ 1.5 x 10 37 erg s -1 , will hide the thermal emission from the cooling ashes. 

The values reported for Model B are qualitatively in agreement with those inferred 
from several XRB sources, such as the textbook burster GS 1826-24 [r rec = 5.74 ± 0.13 h, 
a = 41.7 ± 1.6], 4U 1323-62 [r rec = 5.3 h a = 38 ±41, o r 4U 1 608-52 [r rec = 4.14-7.5 h, 



a 



41-54], to quote a few examples from 



Galloway et al. 



(120081 ) . From the nucleosynthesis 



viewpoint, the main nuclear flow is dominated by the rp-process (rapid proton-captures and 
/3 + -decays), the 3a-reaction, and the ap-process (a sequence of (a, p) and (p, 7) reactions). 
The main flow proceeds away from the valley of stability, eventually reaching the proton 
drip line beyond A = 38. In this context, much of the initial H and 4 He is transformed 
into heavier nuclei. However, it is important to stress the presence of some unburned H 
(X = 0.05) and 4 He (Y = 0.11) in the envelope at the end of the fourth burst (see Table 
[TO]) . Indeed, the mean, mass-averaged metallicity at the end of the fourth bursting episode 
has increased from the original (accreted) metallicity of Zj ~ 0.02 to a value of Z ~ 0.84. 
The nuclei in the envelope with the largest mass fractions are 60 Ni (0.29), 64 Zn (0.14), 4 He 
(0.11), 32 S (0.09), and 56 Ni (0.06). A mass-averaged 12 C yield of ~ 0.02 is obtained at the 
end of the fourth burst, not enough to power a super b urst (which requires X .( 12 C) m .i n ~ 0.1 



at the e nvelope base; see 



Cumm ing 



(120041 ): ICooper k Naravanl (120051 ) 



fc Bildstenl (l200lh : 



Strohmaver fc Brqwnl (20021 ); 



Cumming! (120051 ): ICooper & Naravanl ( 120061 ) ). The 



Brown 



corresponding nucleosynthetic endpoint (defined by the heaviest isotope with a mass 
fraction A,- > 10 -9 ) is found around 100 Pd. 
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Table 7. Summary of burst properties for Model A. 



Burst number T peak (GK) r rec (h) L peak (L ) a 



1 


1.10 


5.95 


1.4 x 10 5 


50 


2 


1.25 


4.63 


2.9 x 10 5 


29 


3 


1.21 


5.51 


2.3 x 10 5 


38 


4 


1.21 


5.21 


2.7 x 10 5 


30 



Table 8. Summary of burst properties for Model B. 



Burst number T peafc (GK) r rec (h) L peak (L ) a 



1 


1.06 


5.88 


9.7 x 10 4 


50 


2 


1.25 


4.44 


2.1 x 10 5 


28 


3 


1.14 


5.19 


1.4 x 10 5 


34 


4 


1.12 


4.84 


1.1 x 10 5 


32 
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Table 9. Summary of burst properties for Model C. 



Burst number T peak (GK) r rec (h) L peak (L Q ) a 



1 


1.03 


5.77 


7.7 x 10 4 


52 


2 


1.26 


4.57 


1.9 x 10 5 


27 


3 


1.17 


5.00 


1.2 x 10 5 


33 


4 


1.17 


4.72 


1.5 x 10 5 


29 
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Table 10. Mean mass fractions of the envelope at the end of the fourth burst, for stable 
nuclides and those with half lives over 1 h and mass fractions X, > 1CT 9 . 



Nucleus Model A Model B Model C 



n 


^ 9 
o.z 


X 


i n- 

1U 


-2 


4 Q 


X 


1 rr 

1U 


-2 


4.4 


X 


i n- 

1U 


-2 


.tie 


Q 7 

y. i 


X 


1U 


-2 


1. 1 


X 


1 n- 

1U 


-1 


l.U 


X 


1U 


-1 


12p 


1.0 


X 


1U 


-2 


1 Q 

i.y 


X 


1U 


-2 


1 7 
1. ( 


X 


1U 


-2 


13 p 


l.o 


X 


ru 


-4 


1.0 


X 


1 n- 
ru 


-4 


l.o 


X 


i ri- 


-4 


14 1\T 
IN 


7 8 
l.o 


X 


i n- 
ru 


-4 


o.4 


X 


1 n- 

1U 


-4 


7 

l.D 


X 


1U 


-4 


IN 


Q 1 

y.o 


X 


ru 


-4 


1 n 

l.U 


X 


1U 


-3 


8 9 


X 


1U 


-4 


16 n 


9 1 
Z. 1 


X 


ru 


-4 


9 fi 


X 


ru 


-4 


9 ^ 
Z.O 


X 


1U 


-4 


17 n 


1. 1 


X 


1 n- 
ru 


-6 


8 9 


X 


1 n- 

1U 


-6 


0.0 


X 


i n- 

1U 


-6 


(j 


1 7 
1. 1 


X 


ru 


-5 


7 n 
( .u 


X 


1U 


-6 


Q 8 

y.o 


X 


ru 


-6 


r 


o.o 


X 


ru 


-5 


1 Q 

i.y 


X 


1U 


-5 


9 ^ 
/.0 


X 


1U 


-5 


19rp 

r 


/I 7 
4. 1 


X 


1U 


-5 


4.0 


X 


1U 


-5 


O.O 


X 


1 n- 


-5 


20 Ne 


1.6 


X 


10" 


-4 


4.2 


X 


10" 


-4 


3.7 


X 


io- 


-4 


21 Ne 


2.5 


X 


10" 


-6 


4.3 


X 


10" 


-6 


4.8 


X 


io- 


-6 


22 Ne 


6.6 


X 


io- 


-6 


2.5 


X 


io- 


-5 


1.8 


X 


io- 


-5 


22 Na 


6.0 


X 


io- 


-4 


6.8 


X 


10" 


-4 


9.0 


X 


io- 


-4 


23 Na 


3.8 


X 


io- 


-5 


1.9 


X 


10" 


-4 


1.2 


X 


io- 


-4 


24 Mg 


4.1 


X 


io- 


-4 


1.4 


X 


10" 


-3 


8.9 


X 


io- 


-4 


25 Mg 


5.0 


X 


io- 


-4 


1.6 


X 


io- 


-3 


1.1 


X 


io- 


-3 


26 Mg 


1.1 


X 


io- 


-3 


1.0 


X 


io- 


-3 


1.1 


X 


io- 


-3 




1.3 


X 


io- 


-5 


1.3 


X 


io- 


-4 


6.6 


X 


io- 


-5 


27 A1 


3.9 


X 


io- 


-4 


7.0 


X 


io- 


-4 


7.8 


X 


10- 


-4 
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Table 10— Continued 



Nucleus Model A Model B Model C 



01 


o.o 


v 

A 


i n - 

1 u 


-3 


1 P, 

-L . U 


v 

A 


1 D" 


-2 


o.o 


v 

A 


1 n 

1 u 


-3 


29c: 


4 ^ 


V 
A 


1 0" 


-4 


7 7 


V 
A 


10" 


-4 


^ 3 

O.O 


A 


1 0" 

X u 


-4 


30 Si 


9 3 


v 

A 


10" 


-3 


3 

O. KJ 


v 

A 


10" 

1W 


-3 


^ 

O. KJ 


v 

A 


10" 


-3 


31p 


1 9 


v 

/\ 


10" 


-3 


1 7 


v 

A 


10" 

1U 


-3 


1 4 

X . t: 


v 

A 


10" 


-3 


32c 
o 


8 7 


v 

A 


10" 

-LU 


-2 


Q 


v 

A 


10" 

XU 


-2 


7 3 


v 

A 


1 0" 

XU 


-2 


33 Q 

kJ 


4 3 


/\ 


lu 


-3 


4 6 

t: . U 


V 


10" 

x u 


-3 


4 

t: . U 


y 


10" 

lu 


-3 


34o 
o 


J. . O 


v 

A 


1 0" 


-2 


1 ^ 
X . o 


v 

A 


1 D" 

-LU 


-2 


1 3 

J. . O 


v 

A 


1 1T 

-LU 


-2 


35 CI 

v_ 1 


1 4 


A 


10" 


-2 


fi 3 
u. o 


V 
A 


10" 

1W 


-3 


7 3 


y 

A 


10" 

XU 


-3 


36 A r 

ill 


X . U 


A 


10" 


-2 


7 

1 . U 


v 

A 


10" 

1W 


-3 


7 7 


y 

A 


10" 

XU 


-3 


37 CI 


3 9 

O. £ 


v 

A 


10" 


-6 


9 ^ 


v 

A 


10" 

1W 


-6 


1 Pi 

X . U 


y 

A 


10" 

XU 


-6 


37 Ar 

ill 


1 9 

X . _ 


v 

A 


10" 


-3 


7 4 


v 

A 


10" 

1W 


-4 


7 fi 


y 

A 


10" 

XU 


-4 


38 A r 


Q ^ 


v 

A 


10" 


-3 


9 


v 

A 


10" 
xu 


-3 


u.o 


y 
A 


10" 

XU 


-3 


39 K 






10" 

lu 


-2 


1 1 

X . X 


y 


10" 

x u 


-2 


1 1 

X . X 


y 


10" 

XU 


-2 


40 Ca 


5.7 


X 


10" 


-3 


4.6 


X 


10" 


-3 


4.0 


X 


io- 


-3 


41 K 


1.2 


X 


10" 


-7 


3.5 


X 


10" 


-8 


3.5 


X 


10" 


-8 


41 Ca 


6.0 


X 


io- 


-5 


7.3 


X 


10" 


-5 


6.3 


X 


io- 


-5 


42 Ca 


2.2 


X 


io- 


-3 


1.5 


X 


10" 


-3 


1.4 


X 


10" 


-3 


43 Ca 


1.4 


X 


io- 


-3 


1.3 


X 


10" 


-3 


1.2 


X 


10" 


-3 


43 Sc 


1.3 


X 


io- 


-3 


8.2 


X 


io- 


-4 


9.5 


X 


10" 


-4 


44 Ca 


7.8 


X 


io- 


-5 


3.1 


X 


10" 


-5 


2.2 


X 


10" 


-5 


44 Sc 


8.8 


X 


io- 


-5 


2.7 


X 


io- 


-5 


4.5 


X 


10" 


-5 


44 Ti 


7.4 


X 


10" 


-4 


6.8 


X 


10" 


-4 


5.7 


X 


10" 


-4 
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Table 10— Continued 



Nucleus Model A Model B Model C 



45o„ 


Q fi 


v 


i n - 

1 u 


-5 


1 "\ 

-L . O 


v 

A 


1 IT 


-4 


8 Q 


v 

A 


i ri" 

-LU 


-5 


45 ry- 

X 1 


7 ^ 


v 

A 


1W 


-5 


Q 7 


V 
A 


1(1" 

1U 


-5 


Q 7 


V 
A 


1(1" 
XU 


-5 


46 T : 


3 3 


x 


10" 


-3 




v 


10" 

x u 


-3 


2 1 

_ . JL 


v 


10" 

X w 


-3 


47 Ti 

X 1 


Q 8 


v 

A 




-4 


8 fi 


v 

A 


1(1" 

1U 


-4 


7 1 

1 . X 


v 

A 


ID" 

XU 


-4 


X 1 


o. o 


v 

A 


10" 

-LU 


-6 


o. z, 


v 

A 


10" 

-LU 


-6 


9 1 


v 

A 


1 0" 
xu 


-6 


48y 

V 


3 8 


/\ 




-4 


3 9 


y 


10" 

x u 


-4 


9 6 

Zj . \J 


y 


10" 

X w 


-4 


48 Cr 


1 Q 

J. . V 


v 

A 


1 fl" 


-3 


9 n 


v 

A 


1 fl" 


-3 


1 7 


y 
A 


1 0" 

-LU 


-3 


49 T : 

X 1 


(\ 7 

U. 1 


v 

A 


i fr 

±u 


-7 


7 Q 


v 

A 


1 0" 
±u 


-7 


7 9 


v 

A 


1 fl" 

-LU 


-7 


49y 

V 


1 1 

J. . J. 


A 


10" 

1W 


-3 


1 4 


V 
A 


1(1" 

-LU 


-3 


1 9 

X . 


y 

A 


If)" 

xu 


-3 


50 n 


3 4 


A 


1W 


-3 


9 Q 


v 

A 


1(1" 

xu 


-3 


9 4 


y 

A 


Ifl" 

xu 


-3 


51y 

V 


3 4 


v 

A 


in - 


-5 


3 ^ 


v 

A 


1(1" 

xu 


-5 


9 7 


y 

A 


Ifl" 

xu 


-5 


51 C 


fi 7 


v 

A 


1(V 


-3 


fi 1 


v 

A 


1 D" 

XU 


-3 




y 

A 


Ifl" 

xu 


-3 


52 Cr 


6 Q 

U . U 


/\ 


10" 


-5 


9 9 


y 


10" 

X w 


-5 


6 5 


y 

A 


10" 

lu 


-5 


52 Mn 


2.9 


X 


10" 


-3 


4.8 


X 


10" 


-3 


3.2 


X 


io- 


-3 


52p e 


6.9 


X 


10" 


-3 


1.2 


X 


10" 


-2 


9.3 


X 


io- 


-3 


53 Mn 


7.7 


X 


10" 


-4 


1.1 


X 


10" 


-3 


7.9 


X 


10" 


-4 


54 Fe 


2.0 


X 


io- 


-3 


1.7 


X 


10" 


-3 


1.3 


X 


10" 


-3 


55 Mn 


1.8 


X 


10" 


-7 


1.5 


X 


10" 


-7 


1.2 


X 


io- 


-7 


55 Fe 


1.2 


X 


io- 


-3 


1.1 


X 


io- 


-3 


8.4 


X 


10" 


-4 


55 Co 


4.7 


X 


io- 


-3 


4.1 


X 


io- 


-3 


3.4 


X 


10" 


-3 


56 Fe 


7.2 


X 


io- 


-6 


3.6 


X 


io- 


-6 


3.2 


X 


10" 


-6 


56 Co 


3.8 


X 


10- 


-3 


2.1 


X 


10" 


-3 


1.9 


X 


10" 


-3 
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56 Ni 

IN 1 


1 4 


A 


i n- 


-1 


U.O 


A 


1U 


-2 


u.u 


V 
A 


1 n 


-2 


57 Fp 


4 7 


A. 




-7 


1 Q 


v 

A 


1 rr 


-7 


-L . U 


v 

A 


1fl" 


-7 


V7 U 


4 7 


A 




-4 


9 fi 


V 
A 


1(1" 

1U 


-4 


1 

x . o 


V 
A 


1(1" 

-LU 


-4 


57 Ni 


R ^ 
o.o 


v 

A 


10" 


-3 


O.O 


v 

A 


10" 

1U 


-3 


9 1 


v 

A 


10" 
xu 


-3 


58 Ni 


A 9 


v 

A 


10" 

-LU 


-3 


9 ^ 


v 

A 


10" 

1U 


-3 


1 

X . u 


v 

A 


10" 
xu 


-3 


59 Ni 

IN 1 


o. o 


v 

A 


10" 


-3 


4 1 

4:. J. 


v 

A 


10" 

1U 


-3 


9 7 


v 

A 


10" 
xu 


-3 


IN 1 


9 ft 


A 


1 0" 


-1 


9 Q 


A 


1 0" 


-1 


9 Q 


A 


1 0" 

1U 


-1 


61 Ni 

IN 1 


1 

l.U 


v 

A 


1 0" 


-3 


9 7 


A 


10" 


-3 


l.O 


V 
A 


10" 

1U 


-3 




1 1 

J- . J- 


A 


10" 


-3 


o.o 


V 
A 


10" 
xu 


-3 


9 8 


A 


10" 
xu 


-3 


62 Ni 

IN 1 


4 3 


v 

A 


10" 


-4 


o.o 


v 

A 


10" 

1U 


-4 


3 1 
0. 1 


v 

A 


10" 
xu 


-4 


62 7,n 


7 Q 


v 

A 


10" 
xu 


-4 


1 7 


v 

A 


10" 


-3 


1 1 

X . 1 


v 

A 


10" 
xu 


-3 


63p 


9 4 


v 

A 


10" 


-3 


3 1 

0. 1 


A 


10" 
1U 


-3 


9 3 


A 


1 0" 
1U 


-3 


64 Zu 


1 4 


A 


10" 


-1 


1 4 


A 


10" 
xu 


-1 


1 7 


V 
A 


10" 
xu 


-1 


65 Cu 


7.5 


X 


10" 


-7 


1.5 


X 


10" 


-6 


1.2 


X 


10" 


-6 


65 Zn 


1.2 


X 


10" 


-3 


2.8 


X 


10" 


-3 


2.4 


X 


10" 


-3 


66 Zn 


2.2 


X 


10" 


-4 


4.0 


X 


10" 


-4 


3.3 


X 


10" 


-4 


66 Ga 


4.7 


X 


io- 


-4 


1.0 


X 


10" 


-3 


8.1 


X 


10" 


-4 


66 Ge 


3.3 


X 


io- 


-4 


8.7 


X 


10" 


-4 


8.4 


X 


10" 


-4 


67 Ga 


8.0 


X 


io- 


-4 


1.4 


X 


io- 


-3 


1.3 


X 


10" 


-3 


68 Ge 


4.4 


X 


io- 


-2 


4.7 


X 


10" 


-2 


5.9 


X 


10" 


-2 


69 Ge 


9.5 


X 


io- 


-4 


2.3 


X 


io- 


-3 


2.3 


X 


10" 


-3 


70 Ge 


2.9 


X 


10- 


-4 


8.1 


X 


10- 


-4 


7.5 


X 


10" 


-4 
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71 As 

no 


4 ^ 


v 

A 


i n - 

±u 


-4 


Q Q 


v 

A 


i rr 


-4 


1 1 


v 

A 


1 fl" 

-LU 


-3 


72 c 


1 4 


V 
/\ 


1 fl" 


-2 


X . u 


V 
A 


1(1" 

XU 


-2 


9 9 


V 
A 


1 fl" 

XU 


-2 


73 Se 


4 8 




ID" 


-4 


1 3 

X . O 


v 


10" 

X u 


-3 


1 4 


V 


10" 

XU 


-3 


74 Se 


9 7 


v 

/\ 


.LU 


-4 


u. o 


v 

A 


1(1" 

XU 


-4 


7 3 


v 

A 


1(1" 

XU 


-4 


75 Rr 




v 

A 


10" 

-LU 


-4 


7 

U. I 


v 

A 


10" 

XU 


-4 


7 7 
i . i 


v 

A 


1 (V 

-LU 


-4 


76 Kr 

XVI 


4 Q 


A 


XU 


-3 


U. X 


v 

A 


10" 

XU 


-3 


Q 9 


V 
A 


1 fl" 

XU 


-3 


77 TO 


o. u 


v 

A 


10" 

XU 


-4 


7 8 


v 

A 


10" 

XU 


-4 


1 
x . u 


v 

A 


10" 

XU 


-3 


78 Kr 


1 Q 


v 

A 


10" 

XU 


-4 


^ 4 


v 

A 


10" 

XU 


-4 


7 9 


v 

A 


10" 

XU 


-4 


79 Kr 

IVI 


1 ^ 


A 


10" 

XU 


-4 


o.o 


V 
A 


10" 

XU 


-4 


4 8 


y 

A 


10" 

XU 


-4 


80c 


J. . O 


A 


1(1" 

XU 


-3 


9 ^ 


v 

A 


1(1" 

XU 


-3 


4 1 

t: . X 


V 
A 


1(1" 

XU 


-3 


XVU 


9 9 

_ . _ 


v 

A 


XU 


-4 


o.o 


v 

A 


1(1" 

XU 


-4 


8 ^ 


v 

A 


1(1" 

XU 


-4 


82 c r 


9 


v 

A 


10" 

XU 


-4 


7 ^ 
I . o 


v 

A 


10" 

XU 


-4 


i i 

X . X 


v 

A 


10" 

XU 


-3 


83Q r 

kJl 


1 7 




XU 


-4 


4 3 

t: . O 


y 


10" 

X u 


-4 


7 6 

1 . U 


y 


10" 

XU 


-4 


84 Sr 


6.8 


X 


10" 


-5 


3.8 


X 


io- 


-4 


4.6 


X 


10" 


-4 


85y 


1.9 


X 


10' 


-4 


4.6 


X 


io- 


-4 


6.9 


X 


10" 


-4 


86 Zr 


2.1 


X 


io- 


-4 


6.3 


X 


io- 


-4 


9.8 


X 


10" 


-4 


87 Zr 


1.6 


X 


io- 


-4 


5.4 


X 


10" 


-4 


9.4 


X 


10" 


-4 


88 Zr 


2.5 


X 


io- 


-4 


2.8 


X 


10" 


-4 


7.3 


X 


10" 


-4 


89 Nb 


6.5 


X 


io- 


-4 


6.1 


X 


10" 


-4 


2.4 


X 


10" 


-3 


90 Mo 


1.2 


X 


io- 


-4 


2.7 


X 


10" 


-4 


1.0 


X 


10" 


-3 


91 Nb 


3.7 


X 


io- 


-5 


1.5 


X 


io- 


-4 


5.8 


X 


10" 


-4 


92 Mo 


9.2 


X 


10- 


-6 


7.1 


X 


10" 


-5 


2.6 


X 


10" 


-4 
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When the lower limit for the 15 0(a,7) 19 Ne reaction rate is adopted (Model A), peak 
temperatures and luminosities of T pea k ~ 1.1 — 1.3 GK and L pea ^ ~ (1 — 3) x 10 5 L & are 
found (see Table [7]). Recurrence times between bursts around r rec ~ 4.5-6 h and ratios 
between persistent and burst luminosities, a ~ 29 — 50, have also been obtained. All in 
all, these values are very similar to those found with Model B. But more important, the 
model exhibits quasi-pe riodic bursting episo des in sharp contrast with the results reported 



in a similar analysis by 



Fisker et al. 



(120061 ) . In that work, the adopted lower limit on 
the reaction rate, which happens to lie within a factor of two from the lower limit found 
in the present Monte Carlo simulation, led to steady state burning of the accreted layer 
without bursting. With respect to nucleosynthesis, similar amounts of unburned H (0.03, 
by mass) and 4 He (0.1) are found at the end of the fourth burst (see Table ITU]) . The mean, 
mass-averaged metallicity at the end of the fourth bursting episode reached a value of 
Z ~ 0.87 in this model. The most abundant nuclides in the envelope in this model are 60 Ni 
(0.26), 64 Zn (0.14), 56 Ni (0.14), 4 He (0.1), and 32 S (0.09). A mass-averaged 12 C yield of 
~ 0.02 is also obtained at the end of the fourth burst, while the nucleosynthetic endpoint 
corresponds to 101 Pd. 

Finally, when the upper limit for the 15 0(a,7) rate is adopted (Model C), peak 
temperatures and luminosities of T pea ^ ~ 1 — 1.3 GK and L pea ^ ~ (0.8 — 2) x 10 5 L Q , 
are found (with the lower value obtained for the first burst, interpreted as driven by the 
initial conditions; see Table [9]). Recurrence times between bursts around r rec ~ 4.5-6 h and 
ratios between persistent and burst luminosities, a ~ 27 — 52, have also been obtained. 
In summary, the results are very similar to those reported above for Model B. As for 
the associated nucleosynthesis, similar amounts of unburned H (0.04, by mass) and 4 He 
(0.1) are found at the end of the fourth burst (see Table [TO]) . The mean, mass-averaged 
metallicity at the end of this fourth bursting episode is Z ~ 0.86 for this model. The most 
abundant species in the envelope are 60 Ni (0.29), 64 Zn (0.17), 4 He (0.1), 32 S (0.07), and 
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93 Tc 


1.4 x 


10" 


-5 


7.7 x 10" 


-5 


2.6 


x 


10" 


-4 


94 Tc 


1.7 x 


10" 


-5 


6.9 x 10" 


-5 


3.0 


x 


10" 


-4 


95 Ru 


7.7 x 


10" 


-7 


8.4 x 10- 


-6 


5.3 


x 


10" 


-5 


96 Ru 


1.7 x 


10" 


-7 


1.5 x 10- 


-6 


7.6 


x 


10" 


-6 


97 Ru 


1.7 x 


10" 


-7 


6.6 x 10" 


-7 


5.8 


x 


10" 


-6 


98 Ru 


5.2 x 


io- 


-8 


1.5 x 10- 


-7 


2.6 


x 


10" 


-6 


"Rh 


1.3 x 


io- 


-8 


2.2 x 10" 


-8 


8.5 


x 


10" 


-7 


ioo pd 


2.4 x 


10" 


-9 


3.3 x 10" 


-9 


2.1 


x 


10" 


-7 


ioi pd 


1.4 x 


10" 


-9 






9.7 


x 


10" 


-8 


102p d 












1.6 


x 


10" 


-8 


103 Ag 












1.5 


x 


10" 


-9 
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56 Ni (0.06). A mass-averaged 12 C yield of ~ 0.02 is also obtained at the end of the fourth 
burst, while the nucleosynthetic endpoint corresponds now to 103 Ag, slightly above the 
endpoints reported for Models A & B. Indeed, as shown in Table [TUj there is a trend when 
comparing the nucleosynthesis associated with Models A, B, & C: a small increase in the 
final abundances of intermediate-mass and heavy elements, likely driven by the specific 
choice of 15 0(a,7) 19 Ne rate, since a larger adopted value for the rate will translate into a 
more effective breakout from the CNO cycles. 



Conclusion 



In summary, we have evaluated the available experimental data on the states in 19 Ne 
that serve as resonances in the 15 0(a,7) 19 Ne reaction in Type I x-ray bursts. Using the 
likelihood distributions inferred from a decay branching ratio and lifetime measurements, 
we calculated the thermally averaged 15 0(a,7) 19 Ne reaction rate via the Monte Carlo 
method, determining the median and 99.73% confidence level upper and lower limits on the 
reaction rate. At this confidence level, the uncertainty in the rate is very large, exceeding a 
factor of 1000 at some temperatures. Given this large unce rtainty, the central value o f the 



rate, which does not differ from the recommended rates of 



Fisker et al. 



Caughlan fc Fowlerl (119881 ) or 



(120061 ) by more than a factor of three at these temperatures, is not particularly 



meaningful and the whole range of reaction rates allowed by experiment must be considered. 

Nevertheless, with the exception of the weak nucleosynthetic trend outlined above, 
we conclude that variation of the 15 0(a,7) 19 Ne reaction rate within the present 99.73% 
confidence level allowed range has no substantial effect on either the burst properties or the 
accompanying nuc leosynthesis in our models of Type I x-ray bursts. In striking contrast to 



Fisker et al. 



( 120061 ). we do not find that a small value of the 0(a,7) Ne reaction rate at 



the lower limit allowed by experiment leads to steady state burning but rather to marginally 
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Fig. 9. — Peak luminosities attained in the first four x-ray bursts and the time they were 
reached in the model calculations using three different 15 0(a,7) 19 Ne reaction rates. 
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more energetic bursts than are predicted with the recommended value of the reaction rate. 
Further astrophysical model calculations are required to resolve the apparent discrepancy 
between different sensitivity studies concerning the importance of the 15 0(o;,7) 19 Ne reaction 
in models of Type I x-ray bursts. 
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